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Abstract 

We investigate a generalized non-linear 0(3) a-model in three space dimen- 
sions where the fields are maps from M 3 U{oo} to S 2 . Such maps are classified 
by a homotopy invariant called the Hopf number which takes integer values. 
The model exhibits soliton solutions of closed vortex type which have a lower 
topological bound on their energies. We numerically compute the fields for 
topological charge 1 and 2 and discuss their shapes and binding energies. The 
effect of an additional potential term is considered and an approximation is 

given for the spectrum of slowly rotating solitons. 
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I. INTRODUCTION 



The non-linear 0(3) cr-model in (3+1) dimensional space-time is a scalar field theory 
whose target space is S 2 . The static fields are maps R 3 U {00} 1— > S 2 and can be classified 
by a homotopy invariant which is called the Hopf number. Such a model in three space 
dimensions must include higher order terms in the field gradient in order to allow non- 
singular, topologically non-trivial, static solutions. The corresponding "sigma model with a 
Skyrme term" was proposed long ago by L.D. Faddeev PJ. For this model the Hopf number 
provides a lower topological bound on the energy || . 

Early studies on "Hopfions" (soliton solutions of Hopf number unity) in classical field 
theory, including estimates for their size and mass, were carried out by de Vega ||. Subse- 
quently it was suggested to employ them in an effective chiral theory describing low-energy 
hadron dynamics; in that respect they are similar to Skyrmions Q. It was later shown by 
Kundu and Rybakov [|J that Hopfions in the 0(3) a- model are of closed vortex type. 

Models with non-zero Hopf number have also been investigated in condensed matter 
physics for the description of three-dimensional ferromagnets and superfluid 3 He |||7|]. These 
are effective theories of Ginzburg-Landau type where the fields are interpreted as physical 
order parameters. However, a field configuration which is a solution of the full equations of 
motion has not been found for any of the mentioned theories. 

In this paper we mainly study classical static Hopfions. Our model is defined in section [H 
where also an ansatz of azimuthal symmetry is introduced which is later used for numerical 



computations. In section |T| we present our numerical results which are minima of the energy 
functional for Hopf number one and two. We discuss their shapes and binding energies as well 
as their relation to (2+1) dimensional solitons. Our model has a self-interaction coupling pa- 
rameter and we study the dependence of the energy on this coupling. In addition, the effect 
of a symmetry breaking potential term is described. In section |Iy| we give a simple approxi- 
mation for the excitation spectrum of a Hopfion slowly rotating around its axis of symmetry. 
We conclude with section [V] where we also remark on possible further investigations. 
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II. HOPF MAPS AND TOROIDAL ANSATZ 



We are almost exclusively interested in static solutions and therefore define our model 
by the following energy functional on R 3 

E s tat [4>] = A / d-K ~ {d^f + ^ (9,0 x <9» 2 + f (<9,</>) 2 (d^) 2 . (1) 

For ^2 = this is equivalent to the static energy of the Faddeev-Skyrme model |j]J2|. The 
field is a three component vector in iso-space, subject to the constraint 2 = 1. The cross 
product is taken in internal space and the coordinate indices i,j run from 1 to 3. 

For gi = g 2 = minima of E (eq. (HD) are harmonic maps from R 3 to S 2 . As shown in 
ref. |§, all non-constant harmonic maps are orthogonal projections R 3 i— > R 2 , followed by a 
harmonic map R 2 i— > S 2 and therefore have infinite energy. 

Consistently, simple scaling arguments along the line of the Hobart-Derrick theorem || 
show that the fourth order terms in the energy functional are required to stabilize the soliton 
against shrinkage. We include here the most general combination of global 0(3)-invariant 
fourth order terms. 

The parameter A is a constant of dimension energy/length and determines the models 
energy unit. The couplings g\ and g-i are of dimension (length) 2 . The ratio g\jg-i is the 
only physically relevant coupling since an overall scaling of g\ and g<i can be absorbed by a 
rescaling of length and energy units. Using (di<f> x dj(p) 2 = {di<f)) 2 {dj<fi) 2 — (di<p ■ dj<fi) 2 and 
the inequality 

2 &<t> ■ d^f > E (<^) 2 ^ E ^ ■ d ^ 2 > ( 2 ) 

ij ij ij 

one sees that the allowed ranges for the coupling constants are g 2 > and g\ > —2g 2 . 

For finite energy solutions one requires <fi — > n as |r| — >• oo, where n is a constant unit 
vector. Thus R 3 can be one-point compactified to S 3 and the fields (ft are maps 

tf>: S 3 ^ S 2 . (3) 
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Because n 3 (S 2 ) = Z, every (j> falls into a class of topologically equivalent maps, where each 
class is characterized by an integer: the Hopf number H. 

Although it is not a simple "winding number", H has an elementary geometric inter- 
pretation. The pre-image of every point of the target space S 2 is isomorphic to a circle. 
All those circles are interlinked with each other in the sense that any circle intersects the 
disc spanned by any other one. The Hopf number just equals the multiplicity by which two 
arbitrary circles are linked. 

H also has a differential geometric representation JlO| : If / is a generator of the de-Rham 



cohomology H^ R (S 2 ), its pullback F under <fi is necessarily exact since Hl, R (S 3 ) = 0. Hence 
a 1-form A with F = dA exist and H ~ / A A F. 

In coordinate language, the dual of F is Bi = e^k 4> • dj<p x dk4> and 

H = -— \, [ g?x B • A . (4) 



(8n) 2 Jr 

It was proved in [[| that the energy eq. (|1|) has a lower topological bound in terms of H . 
For 9i > it is given by 

E stat >AkH 3 /\ (5) 

where k = y/2^(2n) 2 3 3 ^ §. 

The variational equations resulting from eq. ([!]) are coupled non-linear partial differential 
equations. It would be useful to find a parametrization of 4> which carries non-zero Hopf 
charge and allows the equations to be reduced to ordinary differential equations. There have 
been two proposals for such fields in the literature. One of them uses spherical coordinates 
and is a composition of the standard Hopf map and a map S 2 i— > S 2 for which a hedgehog 
ansatz is employed Alternatively, a closed vortex ansatz in toroidal coordinates was 



suggested P,^ JI^JI3| . However, as shown in |TJ]], even for g 2 = none of these proposals 
allows a consistent separation of variables in the variational equations derived from eq. ([[]) . 

At this point it is instructive to look at the symmetries of the field. It was shown in 
ref. that the maximal subgroup of 0(3)x ® 0(3)/ under which fields with non-vanishing 
Hopf number can be invariant is 



G = diag[O(2)x®0(2)j]. (6) 

Here 0(2)x and 0(2)/ denote rotations about a fixed axis in space and iso-space respectively. 
We choose the z- and 03-axis as the axes of symmetry. According to the Coleman-Palais 
theorem we expect to find the minimal energy solution in the class of G-invariant config- 



urations ||17|| . Therefore we use the most general G-invariant ansatz, written in terms of 
two functions w(^i,^ 2 ) and ^(^1,^2)- They depend on coordinates £1 and £ 2 which form an 
orthogonal coordinate system together with a, the angle around the z-axis: 

0i + ifa = V^W^^ Na+V ^ 2)) , 03 = &)• (7) 

We have checked the consistency of this ansatz with the variational equations derived from 
eq. (0). The components 0i and 02 have to vanish along the z-axis for the field to be 
well-defined. This is realized by setting 0(0, 0, z) = n = (0, 0, 1), which also defines the 
vacuum state of the theory. In order to describe a non-trivial map, has to be surjective. 
Hence there is at least one point r with 0(r o ) = — n. Under the action of G, r represents 
a full circle around the z-axis. We fix our coordinate system such that this circle lies in 
the xy-pl&ne and define a = |r |. On every trajectory from the circle to the z-axis or 
infinity, w(£i,£2) runs at least once from —1 to 1. Therefore the surfaces of constant w are 
homeomorphic to tori. 

This structure prompts us to choose toroidal coordinates (77, /3, at), related to cylindrical 
coordinates (r, z, a) as follows 

asmhn a sin/? 
r = '-, z= ^, (8) 

T T 

where r = cosh rj — cos (3. Surfaces of constant rj describe tori about the z-axis, while each 
of these tori is parametrized by the two angles (/3,a). The two cases rj = and rj = 00 
correspond to degenerated tori, rj = being the z-axis and rj = 00 the circle of radius a in 
the xy-pl&ne. 

The function w(r],/3) is subject to the boundary conditions w(0, (3) = l,u>(oo,/3) = —1 
and is periodic in (3. v(r],(3) is an angle around 03 and can include windings around (3. 



Therefore we set v(t),/3) = MP + vo(r],P) where vq(.,P) : S* 1 i— > S 1 is homotopic to the 
constant map. Since v is ill-defined for w = ±1, it is not restricted by any boundary 
condition at 77 = 0,oo. 

The "potential" A and the "field strength" B for this ansatz are given by 

A a = 2 T N(w-1), A p = 2-(M + v Q )(w + l) 1 A v = 2-v' (w + 1) , 
a sinn 77 a a 

2 2 2 

B a = 2 T -{w\M + v G )-v' Q w) ) Bp = -2 T A^', B„ = 2 2 T , iVw, (9) 
a cr sinn 77 er sum 77 

where the dot and prime denote derivatives with respect to (3 and 77 respectively. Note that 
the field A is well defined on all of M 3 . The gauge has been chosen such that A a vanishes 
for rj = (where the coordinate a is undefined) and analogously Ap vanishes for 77 = 00. 

Eq. (|) then gives H = N M in accordance with the linking number definition given 
above. The energy eq. (|l|) of ansatz eq. (0) is given by 

E[w(r], P), v(r), P), a] = ttA [ d V dp + (1 - w") ( (Vv) 2 + ^ 



t 3 1 1 — w 2 \ a 2 sinh 77 

l \ a smn 77 



92 
+ 4 



+ (!-<//) (Vt;) 2 + 



(10) 



1 — w z \ a 2 smh 77 , 

In toroidal coordinates the gradient includes a factor a -1 . Hence the term quadratic in the 
gradients is proportional to a while the quartic terms are inverse proportional to it. For 
soliton solutions, the energy functional has to be varied with respect to w, v and a. 

III. NUMERICAL RESULTS 

The variational equations for eq. (|T0| ) are highly nonlinear coupled PDEs and numerically 
hard to tackle. Therefore we solved the problem by a minimization of the energy functional 
which was discretized on an (77, 0) grid. The search for the minimum in a high-dimensional 



space is feasible using the NETLIB routine t>e08 with an algorithm described in flSfl. This 



method is applicable if the objective function is a sum /(x) = /«( x ) °f simpler functions 



fi, each of which is non-constant only for a few components of the (multi-dimensional) 
vector x. Thus the Hessian matrix is very sparse and can be updated locally This saves 
a considerable amount of memory and time compared to a more naive implementation of a 
conjugate gradient search. 

We obtain field configurations as displayed in Fig. [T|(a) where the Hopf number equals 1. 
In this plot the field (p is viewed from above the north pole of target S 2 . Iso- vectors in the 
northern hemisphere terminate in a cross, those in the southern hemisphere in a dot. The 
toroidal structure of the fields is clearly visible. Also note that the fields in the southern 
hemisphere span a torus indeed. 

There is an interesting interpretation of such configurations in terms of the 0(3) cr-model 
in (2+1) dimensions, the solutions of which we call (anti-) baby Skyrmions. The fields in 
the positive and negative x-halfplane of Fig. [I] are baby Skyrmions and anti-baby Skyrmions 
respectively. This can be understood in the following way. Wilczek and Zee [|1(J] show 



that a (2+l)-dimensional configuration of Hopf number one can be produced by creating 
a baby Skyrmion/anti-baby Skyrmion pair from the vacuum, rotating the (anti-) Skyrmion 
adiabatically by 2ir and then annihilating the pair. In our model time corresponds to the 
third space dimension, hence Fig. [TJ(a) displays a "snapshot" at the time when the anti-baby 
Skyrmion is rotated by ir. Baby Skyrmions are classified by a homotopy invariant Q G Z 
due to 7T2{S 2 ) = Z. The analytic expression for Q is given by 

Q = _L I dx0-0i0 x <9 2 0, (11) 



An _ 

where 1 and 2 denote cartesian coordinates in M 2 . The topological charge density is half the 
a-component of B. The integral over the whole plane vanishes because the contributions for 
negative and for positive x exactly cancel. However, if integrated over the positive half-plane 
only, eq. flTT|) yields the baby Skyrmion number for ansatz (|7|): 

i r 2w r°° a 2 

Q = 7T d P / d V -B a = M, (12) 

o" 71 " -'0 JO T 



where we use B a of eq. 



Next we turn to Hopfions of topological charge two. For parametrisation eq. (0) there 
are two ways of creating a Hopfion with H = 2, namely by setting either N or M to 2. 
Both cases correspond to two Hopfions sitting on top of each other. In order to determine 
which configuration represents the true ground state we computed their energies and found 
that the configuration with N = 2, M = 1 yields the lower energy for all couplings. The 
interpretation of the H = 2 solutions in terms of a (2+l)-dimensional soliton/anti-soliton 
pair is equivalent to the one given above for the 1-Hopfion. Because the multiplicity of 
the azimuthal rotation is N = 2 for the 2-Hopfion, the anti-baby Skyrmion in the negative 
x-halfplane (see Fig. 1(b)) has a relative angle of tt compared to the anti-baby Skyrmion of 
Fig. 1(a). 

It is instructive to investigate how the inclusion of a potential term V[<f>] alters the 
configuration. Its energy can be lowered by rescaling x — > Ax, (A — > 0) under which 
V — > A 3 V^. This means that the potential term induces a "shrinkage" of the configuration in 
the sense that the favoured position of the fields is closer to their vacuum value. This effect 
is counter-balanced by the higher order derivatives in the energy functional eq. ([I]). 

Any potential explicitly breaks the models global 0(3) symmetry because 0(3) acts 
transitively on the target space. We chose V = m 2 J cfx(l — n • </>), where the parameter 
m is of dimension (length) -1 and, in a quantum version of the theory, becomes the mass 
of the elementary excitations. The minimum energy solution for m = 4 can be seen in 
Fig. 1(c). The tube-like region where the field is in the southern hemisphere has clearly 
shrunk. Adding a linear potential term also means that the fields fall off exponentially at 
large distances. The reason is that the equations of motion become in the asymptotic limit 
those of the massive Klein-Gordon equation. 

The fields of minimal energy correspond, via eq. ([I]), to energy distributions which are 
displayed in Fig. 2. Despite the toroidal structure of the fields, we find that the energy 
for the Hopfion of H — 1 is lump-shaped, see Fig. 2(a). Although unexpected, this is not 
entirely unlikely, because the field changes far more rapidly within the disc |r| < a than 
outside it. Hence the gradient energy can be concentrated in the vicinity of the origin. 
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If the potential term becomes very large compared to the gradient terms one expects 
the energy to become more localized around the filament where the fields are far away from 
the vacuum. We observe this transition to a toroidal energy distribution at m ~ 4 for 
9i = ^,92 = 0. 

The energy distribution of the 2-Hopfion is of toroidal shape (for all m), as shown in 
Fig. 2(b). It is a common feature in many soliton theories that solutions of topological charge 
two are tori, notably Skyrmions, baby Skyrmions and magnetic monopoles. It is interesting 
to ask whether the 2-Hopfion is in a stable state or likely to decay into two Hopfions of 
charge one. As an estimate for short range interactions one can compare the energy per 
Hopfion for the solution of H = 1 and H = 2 and conclude from the sign of the energy 
gap whether there is a repulsive or attractive channel. Our results are plotted in Fig. 3(a), 
which also shows the topological bound eq. @. For a pure Skyrme coupling we obtain 
energies of 197A and 2 * 158A for the 1-Hopfion and 2-Hopfion respectively. Moreover, it 
turns out that for all couplings the 2-Hopfion has a lower energy per topological unit than 
the 1-Hopfion. This indicates that there is a range where the forces are attractive and that 
the 2-Hopfion can be stable at least under small perturbations. Of course, there can be a 
range in which the forces are repulsive, however, an investigation of such interactions would 
require a full (3+l)-dimensional simulation which is beyond our present means. Also note 
that the gap between the energies per Hopfion is largest when the fourth order terms are 
purely the Skyrme term. On the other hand, for g% — > —2g 2 , (i.e. g — > 1) the energy of the 
quartic terms tends to zero. Hence the energy of the soliton vanishes as a consequence of 
the above mentioned Hobart-Derrick theorem. 

IV. SPINNING HOPFIONS 

Finally, we study the effect of a slow rotation around the axis of symmetry. For this 
we use a Lorentz-invariant extension of our model into (3+1) dimensional space-time. The 
energy of the rotating Hopfion E = E rot + E stat , where E stat is the static energy given by 
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eq. (|l[) and E rot is the rotation energy functional: 

E rot [(f)] = A [ d-K I (d t <f>) 2 + f (d t <f> x d l( j>) 2 + f (<%0) 2 («9,</>) 2 + O ((^0) 4 ) . (13) 

In the spirit of a moduli space approximation we assume that the configuration does not alter 
its shape due to the rotation ("rigid rotor"), i.e. it is given at any time by a static solution 
(see Jl7|] for a review on similar treatment of the Skyrmion). We impose time dependence 
on the azimuthal angle by a — > a + j^t with constant velocity to. E rot leads to a term in the 
energy that is proportional to uo 2 

E = E stat + J -u 2 , (14) 



1-w 2 ) 



where terms 0(u; 4 ) are neglected. J is the moment of inertia and, using eq. fl7|), given by 

(15) 

J can be measured explicitly on the individual solution. We plotted the values for H = 1 
and H = 2 in Fig. 3(b). The moment of inertia per Hopfion is always larger for the 
H = 1 solution, with an increasing gap for decreasing g. This should be compared with the 
dependence of E stat on g. 

The functional E stat (eq. ([]])) is invariant under a- rotations while the fields of ansatz (0) 
are clearly not. Therefore, upon quantization, the coordinate a describes a zero- mode and 
requires treatment as a collective coordinate. This is similar to the problem of the rotating 
radially symmetric Skyrmion. In analogy to the Skyrme model we therefore USG, clS db first 
approximation, the spectrum obtained by a straightforward quantization. The canonical 
momentum is I = i-4-, (h = 1) and the rotational energy E rot = —l 2 /2J. It is then trivial to 

2 

solve the eigenvalue problem E rot i/j = Xi/j, which gives A n = |j. 

V. CONCLUSIONS 



We have studied topological solitons in a generalized non-linear 0(3) cr-model in three 
space dimensions. Physically one may think of them as a model for hadronic matter or 
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topological defects in a condensed matter system. By using a general ansatz for the fields 
we obtained explicit numerical solutions for soliton number one and two. Unexpectedly, the 
energy of the 1-Hopfion is distributed as a lump. We also observed that two solitons sitting 
on top of each other attract, thus indicating a stable configuration. 

There are several interesting questions which remain unanswered. In particular, the 
stability of Hopfions of higher topological charge deserves some scrutiny. It is worthwhile 
asking how multi-solitons which sit on top of each other, or at least are very close, behave 
under induced perturbations. In analogy to planar 0(3) a-models there might be several 
decay channels into less symmetric configurations |L8] . 



At the opposite end of the scale, it would be instructive to look in greater detail at 
the interaction potential of two or more well-separated Hopfions. This is also interesting in 
comparison to the well-studied dynamics of Skyrmions and monopoles. Clearly, a first step 
in such an investigation would be to determine the asymptotic fields of the Hopf soliton. It 
seems obvious that inter-soliton forces will depend on the orientation of the Hopfions. 

The complete description of Hopfion dynamics would require a huge numerical effort 
which can, however, possibly be reduced by an appropriate approximation scheme. For 
Bogomol'nyi solitons, the low-energy behaviour can be approximated via the truncation of 
the dynamics to the moduli space. Although our numerical results show that Hopfions are 
not of Bogomol'nyi type, given that the static forces between them are weak, there is a 
chance that their dynamics can be described by some kind of moduli space approximation, 
in analogy to Skyrmions (which are also not of Bogomol'nyi type). 

Finally, it seems worth to study spinning Hopfions in a more sophisticated way. This 
should include an assessment of the back reaction of the rotation on the matter fields. From 
this one expects a non-trivial shift of the energy levels in the rotation spectrum and possibly 
radiation of excessive energy. 
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FIGURES 

FIG. 1. (a) Field configuration in the xz-plane for H = l,gi = 0.4, g 2 = 0.4. (b) Field 
configuration of H = 2,g\ = 0.4, g<i = 0.4. (c) Field configuration with potential term, 
H = l,gi = 1,(72 = 0,m = 4. The field is projected into the </>i02-plane. A cross indicates 
4>s > 0, a dot ^3 < 0. Therefore the vacuum state is denoted by a cross only. 

FIG. 2. (a) Energy density e (arbitrary units) for H = l,g\ = 0.4, #2 = in cylindrical 
coordinates r, z. (b) Energy density e for H = 2, g\ = 0.4, gi = 0.8 over r, z. 

FIG. 3. (a) Dependence of the energy E s t a t per Hopfion on the quartic couplings. They are 
parametrized as g\ = 1 — 3g,g2 = 9- Hence g = corresponds to pure (di4> x 3j0) 2 coupling, 
^ = 1/3 to pure (di4>) 2 (dj4>) 2 coupling and g = 1 to the case (di<j)) 2 (dj<j)) 2 — 2 x <9j0) 2 . The 
energy is given in units of A. The topological bounds for pure Skyrme coupling are also displayed, 
(b) Dependence of the moment of inertia J per Hopfion on the coupling g. 
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Fig 1(a) 
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